Novel scale-free small-world networks from Koch curves 
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The class of Koch fractals is one of the most interesting families of fractals, and the study of 
complex networks is a central issue in the scientific community. In this paper, inspired by the 
famous Koch fractals, we propose a mapping technique converting Koch fractals into a family 
of deterministic networks, called Koch networks. This novel class of networks incorporates some 
key properties characterizing a majority of real-life networked systems — a power-law distribution 
with exponent in the range between 2 and 3, a high clustering coefficient, small average path 
length, and degree correlations. All these features are obtained exactly according to the proposed 
generation algorithm of the networks considered. The network representation approach could be 
used to investigate the complexity of some real-world systems from the perspective of complex 
networks, and our networks may describe various real-life networked systems. 

PACS numbers: 89.75.Hc, 05.10.-a, 89.75.Fb, 61.43.Hv 



I. INTRODUCTION 

The past decade has witnessed a great deal of activity 
devoted to complex networks by the scientific community, 
since many systems in the real world can be described and 
characterized by complex networks [l], 0, H, 0] • Prompted 
by the computerization of data acquisition and the in- 
creased computing power of computers, researchers have 
done a lot of empirical studies on diverse real networked 
systems, unveiling the presence of some generic proper- 
ties of various natural and manmade networks: power- 
law degree distribution P(k) ~ fc -7 with characteristic 
exponent 7 in the range between 2 and 3 [5[ , small- world 
effect including large clustering coefficient and small aver- 
age path length ( APL) @ , and degree correlations 0, § • 

The empirical studies have inspired researchers to con- 
struct network models with the aim to reproduce or ex- 
plain the striking common features of real-life systems p], 
2]. In addition to the seminal Watts-Strogatz's (WS) 
small- world network model [6| and Barabasi- Albert's 
(BA) scale-free network model [5| , a considerable number 
of models and mechanisms have been developed to mimic 
real- world systems, including initial attractiveness [9j, 
aging and cost [io| . fitness model fijl, duplication [l2l j. 
weight or traffic driven evolution Il3l. H3 |. geographical 
constraint , accelerating growth [10] , coevolution [ljj , 
visibility graph jl8| . to name but a few. Although sig- 
nificant progress has been made in the field of network 
modeling and has led to a significant improvement in 
our understanding of complex systems, it is still a funda- 
mental task and of current interest to construct models 
mimicking real networks and reproducing their generic 
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properties from different angles. 

In this present article, inspired by the famous class 
of Koch fractals, we propose a family of deterministic 
mathematical networks, called Koch networks, which in- 
tegrates the observed properties of real networks in a sin- 
gle framework. We derive analytically exact scaling laws 
for the degree distribution, the clustering coefficient, the 
average path length, and the degree correlations. The ob- 
tained precise results show that Koch networks have rich 
topological features: they obey power-law degree distri- 
bution with exponent lies between 2 and 3; they have 
a large clustering coefficient, their average path length 
grows logarithmically with the total number of nodes; 
and they may be either disassortative or uncorrelated. 

This work unfolds an alternative perspective in the 
study of complex networks. Instead of searching gen- 
eration mechanisms for real networks, we explore deter- 
ministic mathematical networks that exhibit some typical 
properties of real-world systems. As the classical Koch 
fractals are important for the understanding of geomet- 
rical fractals in real systems |l9j, we believe that Koch 
networks could provide valuable insights into real-world 
systems. 



II. NETWORK CONSTRUCTION 

In order to define the networks, we first introduce a 
classical fractal — Koch curve, which was proposed by 
Helge von Koch [2fJ. The Koch curve, denoted as S±(i) 
after t generations, can be constructed in an recursive 
way. To produce this well-known fractal, we begin with 
an equilateral triangle and let this initial configuration 
be Si(0). In the first generation, we perform the follow- 
ing operations: firstly, we trisect each side of the initial 
equilateral triangle; secondly, on the middle segment of 
each side, we construct new equilateral triangles whose 



FIG. 1: (Color online) The first two generations of the con- 
struction for Koch curve. 



interiors lie external to the region exclosed by the base 
triangle; thirdly, we remove the three middle segments 
of the base triangle, upon which new triangles were es- 
tablished. Thus, we get Si(l). In the second generation, 
for each line segment in Sl(1), repeat above procedure of 
three operations to obtain Si (2). This process is then re- 
peated for successive generations. As t tends to infinite, 
the Koch curve is obtained, and its Hausdorff dimension 
is df = ^ (H- Figured] depicts the structure of Si (2). 

Koch curve can be easily generalized to other dimen- 
sions by introducing a parameter m [2l|. Denoted by 
S m (t) the generalization after t generations, which is con- 
structed as follows [2l| : Start with an equilateral triangle 
as the initial configuration S m (0) . In the first generation, 
we perform the following operations similar to those de- 
scribed in last paragraph: partition each side of the initial 
triangle into 2m + 1 (m is a positive integer) segments, 
which are consecutively numbered 1, 2, • • ■ , 2m, 2m + 1 
from one endpoint of the side to the other; construct 
a new small equilateral triangle on each even-numbered 
segment so that the interiors of the new triangles lie in the 
exterior of the base triangle; remove the segments upon 
which triangles were constructed. In this way we obtain 
S m (l). Analogously, we can get S m (t) from S m (t — 1) by 
repeating recursively the procedure of above three oper- 
ations for each existing line segments in generation t — 1. 
In the infinite t limit, the Hausdorff dimension of the gen- 
eralized Koch curves df = j^m+i) P^l - Figure [H shows 
the the structure of Sa(2). 

The generalized Koch curves can be used as a basis of 
a new class of networks: sides (excluding those deleted) 
of the triangles constructed at arbitrary generations are 
mapped to nodes, which are connected if the correspond- 
ing sides are in contact. For uniformity, the three sides of 
the initial equilateral triangle of S m (0) also correspond 
to three different nodes. We shall call the resultant net- 
works Koch networks. Note that after the birth of each 
side of a triangle constructed at a given generation, al- 



FIG. 2: (Color online) The first two generations of the con- 
struction for generalized Koch curve in the case of m = 2. 




FIG. 3: (Color online) The network derived from Si (2). 



though some segments of it will be removed at subsequent 
steps, we look on its remain segments as a whole and map 
it to only one node. Figures [3] and [4] show two networks 
corresponding to Si (2) and Sa(2), respectively. 

Obviously, Koch networks have an infinite number of 
nodes. But in what follows we shall generally consider the 
network characteristics after a finite number of genera- 
tions in the development of the complete Koch networks. 
From our analytical results, we can quickly obtain the 
characteristics of the complete networks by taking the 
limit of large t. However, the numerical results are nec- 
essarily limited to networks with finite order (number of 
all nodes). 



III. GENERATION ALGORITHM 

According to the construction process of the general- 
ized Koch curves and the proposed method of mapping 
from Koch curves to Koch networks, we can introduce 
with ease an iterative algorithm to create Koch networks, 
denoted by K m . t after t generation evolutions. The al- 
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FIG. 4: (Color online) The network corresponding to S l 2(2). 




and 9m new edges at step t, then one can easily ob- 
tain: L v (t) = 6mL/\(t — 1) = 6m (3m + and 
L e (t) = 9m La (t - 1) = 9m (3m + l)*" 1 , both of which 
hold for arbitrary t > 0. Then, the total number of nodes 
N t and edges E t present at step t is 



t 

E 



L V (U) = 2 (3m + 1)* + 1 



and 



t 

E 

ti=0 



L e (t()= 3 (3ro + l)* 



respectively. Thus, the average degree is 

6(3m+ l) f 



(fe) 



2_E* _ 
N t 2 (3m + !)* + !' 



(1) 



(2) 



(3) 



which is approximately 3 for large i, showing that Koch 
networks are sparse as most real-life networks. 



IV. TOPOLOGICAL PROPERTIES 

Now we study some relevant characteristics of the Koch 
networks K m t , focusing on degree distribution, cluster- 
ing coefficient, average path length, and degree correla- 
tions. 



A. Degree distribution 



FIG. 5: (Color online) Iterative construction method for the 
Koch networks. 



gorithm is as follows: Initially (i = 0), K m ,o consists 
of three nodes forming a triangle. For t > 1, K m j 
is obtained from K m j—\ by adding m groups of nodes 
for each of the three nodes of every existing triangles in 
K m> t-i- Each node group has two nodes. These two new 
nodes and its "mother" node are linked to one another 
shaping a new triangle. In other word, to obtain K m ^ 
from K m .t-i, we replace each of the existing triangles of 
K m j—i by the connected clusters on the right hand of 
Fig. \E\ Figures [3] and [4] illustrate the growing process 
of the networks for two particular cases of m = 1 and 
?7i = 2, respectively. 

Let us compute the order and size (number of all edges) 
of Koch networks K m j. To this end, we first consider 
the total number of triangles L/\(t) that exist at step t. 
By construction (see Fig. O, this quantity increases by 
a factor of 3m + 1, i.e., L&{t) = (3m + 1) L&{t — 1). 
Considering L\{Q) = 1, we have L\{t) — (3m + 1)'. 
Denote L v (t) and L e (t) as the numbers of nodes and 
edges created at step t, respectively. Note that each 
triangle in K m t -i will give birth to 6m new nodes 



Let ki(t) be the degree of a node i at time t. When 
node i enters the network at step ti (ti > 0), it has a 
degree of 2, viz. fcj(tj) = 2. To determine fcj(t), we first 
consider the number of triangles involving node i at step 
t that is denoted by L&(i,t). These triangles will give 
rise to new nodes linked to the node i at step t+1. Then 
at step ti, L&(i, ti) = 1. By construction, for any triangle 
involving node i at a given step, it will lead to m new 
triangles passing by node i at next step. Thus, La(*, t) = 
(m + 1) L&(i,t — 1). Considering the initial condition 
LA{i,ti) = 1, we have LA(i,t) = (m + l) t_ti . On the 
other hand, every triangle passing by node i contains two 
links connected to i, therefore we have ki(t) = 2LA(i,t). 
Then we obtain 



ki{t) =2L A {i,t) = 2(m+ 1) 



t-ti 



(4) 



In this way, at time t the degree of arbitrary node i 
of Koch networks has been computed explicitly. From 
Eq. (|3]), it is easy to see that at each step the degree of 
a node increases m times, i.e., 



h{t) = (m+l)ki(t-l). 



(5) 



Equation (0| shows that the degree spectrum of Koch 
networks is discrete. Thus, we can get the degree dis- 
tribution P(k) of the Koch networks via the cumulative 
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degree distribution 0] given by 

P c Uk) = - V L v (r) = 2X ( 3m + 1 ^ + 1 . (6) 
cumV ' N t f^ t 2x(3m+l) f + l V; 

ln( — ) 

Substituting for ti in this expression using ti = t— ^^pR 
gives 

2x(3m+irx(|)-^ 1 +l 
cumW_ 2x(3m + l)* + l ' li 

In the infinite t limit, we obtain 

ln(3m+l) ln(3m + l) 

Pcum(fc) = 2 x fe EPSTU . (8) 

So the degree distribution follows a power-law form 
P(k) ~ fc~ 7 with the exponent 7 = l+ 1 ^?^ 1 y belonging 
to the interval [2,3]. When m increases from 1 to infi- 
nite, 7 decreases from 3 to 2. It should be stressed that 
the exponent of degree distribution of most real scale- free 
networks also lies in the same range between 2 and 3. 



B. Clustering coefficient 

By definition, the clustering coefficient of a node 
i with degree ki is the ratio between the number of tri- 
angles &i that actually exist among the ki neighbors of 
node i and the maximum possible number of triangles 
involving i, ki(ki — l)/2, namely, Cj = 2e. i /[fc i (fc i — 1)]. 
For Koch networks, we can obtain the exact expression of 
clustering coefficient C(k) for a single node with degree 
k. By construction, for any given node having a degree 
k, there are just & = \ triangles connected this node, see 
also Eq. (HJ). Hence there is a one-to-one corresponding 
relation between the clustering coefficient of a node and 
its degree: for a node of degree k, 



C(k) 



1 



1 



(9) 



which shows a power-law scaling C(k) ~ in the large 
limit of k, in agreement with the behavior observed in a 
variety of real-life systems [22[ • 

After t step growth, the average clustering coefficient 
Ct of the whole network K m ^ tl defined as the mean of 
C i s over all nodes in the network, is given by 



c > = wX 



r=0 



1 



1 



x L v (r) 



(10) 



where the sum runs over all the nodes of all generations 
and G r is the degree of those nodes created at step r, 
which is given by Eq. ([4]). In the limit of large Nt, Equa- 
tion (|10p converges to a nonzero value C, as reported in 
Fig.H For m= 1, 2, and 3, C are 0.82008, 0.88271, and 
0.91316, respectively. As m approaches infinite, C con- 
verges to 1. Thus, C increases with m: when m grows 
from 1 to infinite, C increases form 0.82008 to 1. There- 
fore, for the full range of m, the the average clustering 
coefficient of Koch networks is very high. 
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FIG. 6: Semilogarithmic plot of average clustering coefficient 
Ct versus network order Nt. 




FIG. 7: (Color online) Second construction method of Koch 
networks that highlights self-similarity. The graph after t + 
1 construction steps, K m ,t+i, consists of 2m + 1 copies of 
K m ,t denoted as K^ it (8 = 1, 2, 3, • ■ • , 3m, 3m + 1), which are 
connected to one another as above. 



C. Average path length 

Using a method similar to but different from those in 
literature [23|, [24], [25| , we now study analytically the av- 
erage path length (APL) d t of Koch networks K m , t - It 
follows that 



N t (N t - l)/2 ' 



(11) 



where D t is the total distance between all couples of 
nodes, i.e., 



E 

ieK m , t ,jeK„ 



(12) 



in which rfjj is the shortest distance between node i and 
3- 

Notice that Koch networks have a self-similar struc- 
ture, which allows us to address dt analytically. This 
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self-similar structure is obvious from an equivalent net- 
work construction method: to obtain K m tl one can make 
3m + 1 copies of K mi t-i an d join them in the hubs 
(namely nodes with largest degree) . As shown in Fig. 
network K m j+i may be obtained by the juxtaposition of 
3m + 1 copies of K m . t , which are labeled as t , K^ t , 
• • • , K%%, and K%% +1 , respectively. 

We continue by exhibiting the procedure of the deter- 
mination of the total distance and present the recurrence 
formula, which allows us to obtain Dt+i of the t + 1 gen- 
eration from D t of the t generation. From the obvious 
self-similar structure of Koch networks, it is easy to see 
that the total distance Dt+i satisfies the recursion rela- 
tion 



D t+ i = (3m + 1) D t + Q 



t , 



(13) 



where fit is the sum over all shortest paths whose end- 
points are not in the same t branch. The solution of 
Eq. fHi is 



t-i 



D t = (3m + l) t ~ 1 D l + J2( 3m + r 



f-T-l 



(14) 



T=l 



The paths contributing to f2 t must all go through at least 
one of the three edge nodes (i.e., grey nodes X, Y and 
Z in Fig. [7]) at which the different Kf n t branches are 
connected. The analytical expression for Vl u called the 
length of crossing paths, is found below. 

Let fi"'' 3 be the sum of the length of all shortest paths 
with endpoints in t and t . According to whether 
or not two branches are adjacent, we sort the crossing 
path length f2" ,/3 into two classes: If K^ t and K^ t meet 

at an edge node, ft"' 13 rules out the paths where either 
endpoint is that shared edge node. For example, each 

1 2 

path contributed to fl t ' should not end at node X. If 
t and if^,* do not meet, Slj excludes the paths 
where either endpoint is any edge node. For instance, 
each path contributed to fl 2 ' m+2 should not end at nodes 
X or Y. We can easily compute that the numbers of 
the two types of crossing paths are 3m ^ 3m and 3m 2 , 
respectively. On the other hand, any two crossing paths 
belonging to the same class have identical length. Thus, 
the total sum Qt is given by 



n, 



3m 2 + 3m i 2 



i 2^ ' + Sm IL 



In order to determine fl]' 2 and Sl 2 ' m+2 , we define 



E - 

ieK mit ,i^X 



HX 



(15) 



(16) 



Considering the self-similar network structure, we can 
easily know that at time t + 1, the quantity St+\ evolves 
recursively as 



s t+ i = (m + 1) s t + 2m [s t + (N t - 1)] 
= (3m + 1) s t + 4m (3m + 1)*. 



Using s = 2, we have 

s t = {Ami + 6m + 2) (3m + l)* -1 . 



(18) 



Have obtained St, the next step is to compute the quan- 
tities n]' 2 and Q 2 ' m+2 given by 



n]> 2 = 



E 



^2 (dix+d-jx) 
{Nt-1) d iX + (N t -l) 

2(N t - 1) ]T d lX 

ieifj, t i^x 

2(N t -l)s t , 



3+X 



(19) 



and 



2,m+2 



E di i 



E (* 



X T dXY T UjY , 



ieK^ :t ,ijiX 

2{N t -l)s t + {N t -l) 2 



where dxY = 1 has been used. Substituting Eqs. (|P91 
and ([20]) into Eq. (fl"5|) . we obtain 



(20) 



(17) 



(l t = (9m 2 + 3m)(N t - 1) s t + 3m 2 (N t - l) 2 

= 3m(4m< + 7m + 2)(3m+l) 2t . (21) 

Inserting Eqs. (f2Tj) for fl T into Eq. (Tl4|) . and using D\ — 
53m + 19, we have 

(3m + l)*^ 1 

D t = ± —!- [3m + 5 + (24f + 24m + 4)(3m + if] . 

3 (22) 
Inserting Eq. (J22J) into Eq. (fTTj) . one can obtain the ana- 
lytical expression for dt- 

_ 6m + 10 + 8(6t + 6m+l)(3m+l)* 
' ~ 3(3m+l)[(3m+l) 4 + l] ' ( ' 

which approximates g^f+i in the infinite t, implying that 
the APL shows a logarithmic scaling with network order. 
Therefore, Koch network exhibits a small- world behavior. 
We have checked our analytic result against numerical 
calculations for different m and various t. In all the cases 
we obtain a complete agreement between our theoretical 
formula and the results of numerical investigation, see 
Fig. E 

Recently, it has been suggested that for scale-free net- 
works with order N their average path length d(N) 
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FIG. 8: Average path length dt versus network order Nt on a 
semi-log scale. The solid lines are guides to the eyes. 



is relevant to the exponent 7 [261 |27[ |: when 7 = 3, 
d(N) ~ In TV; when 2 < 7 < 3, d(N) ~ In bxN. However, 
our above exact results for APL of Koch networks show 
that the scaling observed in the literature [H, Ht| is not 
a generic feature of all scale-free networks, at least it is 
not valid for the Koch networks. 



D. Degree correlations 

Degree correlation is a particularly interesting subject 
in the field of network science @, H SI, E3, IM S3 , be- 
cause it can give rise to some interesting network struc- 
ture effects. An interesting quantity related to degree 
correlations is the average degree of the nearest neigh- 
bors for nodes with degree k, denoted as k nn (k), which 
is a function of node degree k [29L [30(| . When k nn (k) in- 
creases with k, it means that nodes have a tendency to 
connect to nodes with a similar or larger degree. In this 
case the network is defined as assortative [JIB]. In con- 
trast, if k nn (k) is decreasing with k, which implies that 
nodes of large degree are likely to have near neighbors 
with small degree, then the network is said to be disas- 
sortative. If correlations are absent, k nn (k) = const. 

We can exactly calculate k nn (k) for Koch networks us- 
ing Eqs. (T?J) and ([5]) to work out how many links are made 
at a particular step to nodes with a particular degree. By 
construction, we have the following expression [31, [33[ 



gree (i.e., t- < i,-) when the node was generated at t%. 
The second sum describes the links made to the current 
smallest degree nodes at each step t\ > ti . The last term 
1 accounts for the link connected to the simultaneously 
emerging node. In order to compute Eq. (|24p . we distin- 
guish two cases according to parameter to: to = 1 and 
to > 2. 

When m = 1, we have 



knn{k) 



L v (U)k(t t ,t) 



f 4 =t 



■ t' i= ti-\ 

J2 mL^kit'iA-Wi^t) 



t'=o 



k nn (k) = t + 2. 



(25) 



Thus, in the case of to = 1, the networks show absence 
of correlations in the full range of t. From Eqs. (|25[) 
and (P) we can easily see that for large t, k nn (k) is ap- 
proximately a logarithmic function of network order Nt, 
namely, k nn (k) ~ lnA^, exhibiting a similar behavior as 
that of the BA model [30l| an d the two-dimensional ran- 
dom Apollonian network [3l| . 

When to > 2, Equation (|24[) is simplified to 



k n n{k) — 



3to + 1) 



TO — 1 

2to 



(to + l) 5 



m - 



1 



3to 
(t - h). 



TO 



(26) 



Thus after the initial step k nn (k) grows linearly with 
time. Writing Eq. in terms of k, it is straightfor- 

ward to obtain 



k n n{k) 



3to+ 1 

TO — 1 



(to + 1) 2 



3m 



- 1 
3 



k 
2 

2 TO 



ln(§ 



TO — 1 



+ 1 ln(m + 1)' 



(27) 



Therefore, k nn (k) is approximately a power law function 
of k with negative exponent, which shows that the net- 
works are disassortative. Note that k nn (k) of the Internet 
exhibits a similar power-lawjiependence on the degree 



with uj = 0.5 m 



CONCLUSIONS 



+ J2 mL v {U)k{uX-l)Wi,t))+l ( 24 ) 

for k = 2 (to+I)* - ** . Here the first sum on the right-hand 
side accounts for the links made to nodes with larger de- 



Networks having a complicated enough structure to 
display nontrivial properties for real-life systems and sta- 
tistical models but simple enough to reveal analytical in- 
sights are few and far between. In this paper, we have 
presented a mapping that converts the Koch curves into 
complex networks, which have many important general 
properties observed in real networks. Our networks are 
characterized by closed-form, exact formulas for various 
properties. The rigorous solutions show that the result- 
ing graphs have a heavy-tailed degree distribution with 
general exponent 7 £ [2,3], a logarithmical average path 
length with network order, large clustering coefficient, 
and degree correlations. Our analytical technique could 
guide and shed light on related studies for deterministic 
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network models by providing a paradigm for computing 
the structural features. 

In addition to the analytical method and the nontriv- 
ial topological characterizations of our proposed model, 
another contribution of our work is mapping the Koch 
fractals to a family of graphs. With the the proposed 
mapping method, a natural bridge between the theory of 
network science and the classic fractals has been built. 
This network representation technique could find appli- 
cation to some real systems making it possible to explore 
the complexity of real-life networked systems in biologi- 
cal and information fields within the framework of com- 
plex network theory. Recently, a similar recipe has been 
adopted for investigating the navigational complexity of 
cities [34(, which has also proven useful to the study of 
polymer physics [35|. Finally, it should be stressed that 
for the special case of m = 1, the network is completely 
uncorrelated, it is thus of potential interest as a null sub- 
strate network to corroborate the dynamical behavior of 



complex systems, since the analytic solution of dynamical 
processes is usually available only for uncorrelated net- 
works [H, [I?], [H, [11] , the generation algorithm of which 
is often more difficult than one would expect a priori [40| . 
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